Raman spectra of BN-nanotubes: Ab-initio and bond-polarizability model calculations 
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We present ab-initio calculations of the non-resonant Raman spectra of zigzag and armchair 
BN nanotubes. In comparison, we implement a generalized bond-polarizability model where the 
parameters are extracted from first-principles calculations of the polarizability tensor of a BN sheet. 
For light-polarization along the tube-axis, the agreement between model and ab-initio spectra is 
almost perfect. For perpendicular polarization, depolarization effects have to be included in the 
model in order to reproduce the ab-initio Raman intensities. 
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Besides its success in the characterization of a large 
range of materials^, Raman spectroscopy has also de- 
veloped into an invaluable tool for the characterization 
of nanotubes. Since the first characterization of (dis- 
ordered) carbon nanotube (CNT) samples^, the tech- 
nique has been refined, including, e.g., polarized Ra- 
man studies of aligned nanotubes^ and isolated tubes^. 
On the theoretical side, non-resonant Raman intensi- 
ties of CNTs have been calculated within the bond- 
polarizability model^. The empirical parameters of this 
model are adapted to fit experimental Raman intensi- 
ties of fullerenes and hydrocarbons. However, the trans- 
ferability of the parameters and the quantitative perfor- 
mance in nanotubes, in particular distinguishing between 
metallic and semiconducting tubes, is still not clear. 

In this communication, we report on the Raman 
spectra of boron nitride nanotubes (BNNTs) 7,8 . Re- 
cently, synthesis of BNNTs in gram quantities has been 
reported 9 . Their characterization through Raman and 
infrared spectroscopy is expected to play an important 
role. However, due to difficulties with the sample purifi- 
cation no experimental data on contamination-free sam- 
ples has been reported. Ab-initio — and empiricaliiii^ 
phonon calculations have determined the position of the 
peaks in the spectra. However, due to missing bond- 
polarizability parameters for BN, the Raman intensi- 
ties have been so far addressed using the model bond- 
polarizability parameters of carbonic. Only the intensi- 
ties of high-frequency modes was presented as it was ar- 
gued that the intensity of low frequency modes are very 
sensitive to the bond-polarizability values^. Here, we 
derive the polarizability parameters for BN sp 2 bonds 
from a single hexagonal BN sheet by calculating the 
polarizability-tensor and its variation under deformation. 
We compare the resulting spectra for BNNTs with full 
ab-initio calculations. We derive conclusions about the 
general applicability of the bond polarizability model for 
semiconducting CNTs. 

In non resonant first-order Raman spectra, peaks ap- 
pear at the frequencies u>„ of the v optical phonon with 
null wave vector with intensities I" which in the Placzek 



approximation^ are 
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Here ej (e s ) is the polarization of the incident (scattered) 
light and n v = [exp(HLJi/ / 'fesT) — with T being the 
temperature. The Raman tensor A" is 
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where is the fcth Cartesian component of atom 7 of 



the vth orthonormal vibrational eigenvector and M 1 
the atomic mass. 
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where £ is the total energy of the unit cell, E is a uniform 
electric field and Uk-y are atomic displacements. This is 
equivalent to the change of the electronic polarizability of 
a unit cell, cvy = Q-Xij (where Q is the unit cell volume 
and Xij t ne electric susceptibility), upon the displace- 
ment Uk<y- The phonon frequencies and eigenvectors^ 
are determined by density functional perturbation the- 
— as implemented in Ref^. For the determination 



oryi 



of the derivative tensor B i J we proceed in two ways: i) 
we calculate it from first principles using the approach of 
Ref. IT5L ii.) we develop a generalized bond-polarizability 
model. 

The basic assumption of the bond polarizability 
modeliiSiii is that the total polarizability can be mod- 
eled in terms of single bond contributions. Each bond is 
assigned a longitudinal polarizability, ai , and a polariza- 
tion perpendicular to the bond, a p . Thus, the polariz- 
ability contribution cqa of a particular bond b is 
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where R is a unit vector along the bond. The sec- 
ond assumption is that the bond polarizabilities only de- 
pend on the bond length R. This allows the calculation 
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TABLE I: Parameters of the bond polarizability model ex- 
tracted from ab-initio calculations (see text). 
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FIG. 1: Unit cell (marked by dashed line) of BN-sheet for the 
calculation of the bond polarizability parameters: a) equilib- 
rium geometry, b) geometry with one bond elongated. 



of the derivative with respect to atomic displacement, 
doqj /du/py, in terms of four parameters ai(R), a p (R), 
a[{R), and a' p (R) (see, e.g. Ref. 0). The use of only 
one perpendicular parameter a p implicitly assumes cylin- 
drical symmetry of the bonds. That can be justified 
in a sp 3 bonding environment. However, in the highly 
anisotropic environment in a sheet of sp 2 bonded carbon 
or BN and the corresponding nanotubes this assumption 
seems hardly justified. In our model we therefore define 
a generalized polarizability with an in-plane (a P i) and 
out-of-plane value (a pa ) of a p . 

With the larger set of parameters, the polarizability 
tensor takes on the more general form 



a i0 = aiRiRj 



(5) 



where S is a unit-vector pointing perpendicular to the 
bond in plane, and T pointing perpendicular to the bond 
out of plane. (In the case of a P i — a po , Eq. (JSJ simplifies 



to Eq. (g} due to the relation S l S j +T i T j 



RiRj.) 



For the derivative tensor (of a single bond), we obtain 



da° 



du 



a't&iRjRk + ai((d k Ri)Rj + MdkRj)) 

+ a' pi SiSjRk + a pi {{d k Si)Sj + Si(d k Sj)) 
+ a'p^fjRk + apoiidkf^fj+fiidhfj)). (6) 

The total derivative tensor is then just the sum over 

all da^j/duk-f of all bonds of the system. The orientation 
of the plane at the position of a particular atom is thereby 
defined by the three nearest neighbor atoms. 

In order to determine the six parameters of our model, 
we perform ab-initio calculations of the polarizability 
tensor ofy of a unit-cell of a single BN-shee1*i2iiS (see 
Fig. ^ a ) ■ The geometry of the system leads to the re- 
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(with the z-axis perpendicular to the sheet). Displacing 
atom 2 in y-direction yields the relation da xx /du2 y 
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FIG. 2: Perpendicular (7±) and longitudinal (711) polarizabil- 
ities per unit length of different BN-nanotubes: ab-initio and 
our generalized bond polarizability model. The influence of 
depolarization can be seen for 7±. 



(3/4) (o{ + a' pl ) + (3/2) (aj + a pl )/R. Finally, by chang- 
ing the geometry of the unit-cell such that one bond is 
elongated while the other two bond-lengths and all the 
bond angles are kept constant (see Fig.^b), we extract 
the derivatives of the bond polarizabilities: a[ — a' xx , 
a' pi = a' yy , and a' po — a' zz . The resulting parameters 
are displayed in Tab. I and compared to the parameters 
we calculated for cubic BN and diamond. The longitudi- 
nal bond polarizability u\ is considerably larger than a p 
which can be intuitively explained as a consequence of 
the "enhanced mobility" of the electrons along the bond. 
For the sheet, the perpendicular polarizabilities clearly 
display different values in the in-plane and out-of-plane 
directions. Without the added flexibility of different pa- 
rameters, the bond-polarizability model would lead to 
inconsistencies in the description of and its deriva- 
tives. In the sheet, an is about twice as large as in cubic 
BN (c-BN) due to the additional contribution of the 7r 
electrons to the longitudinal polarizability. Comparison 
of c-BN with the isoelectronic diamond shows a slightly 
higher polarizability of the C-C bond. 

As a first application of the generalized bond- 
polarizability model, we present in Fig. [21 the polariz- 
ability 7 (per unit length) of different BNNTsS. For the 
polarizability along the tube axis (z-direction) , the model 
(Eq. J5j) agrees almost perfectly with our ab-initio calcu- 
lations. The polarizability is proportional to the number 
of bonds in the unit-cell which is proportional to the tube 
radius. For the perpendicular direction, the model cal- 
culations overestimate the ab-initio values considerably. 
This discrepancy demonstrates the importance of depo- 
larization effects in the perpendicular direction: due to 
the inhomogeneity of the charge distribution in this di- 
rection, an external field induces local fields which coun- 
teract the external field and thereby reduce the overall 
polarizability. The size of this effect can be estimated 
from a simple model. Imagine a dielectric hollow cylin- 
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der of radius R (measured at the mid-point between 
inner and outer wall) and thickness d. The dielectric 
constant in tangential direction, en = (d + i%(3n )/d, is 
different from the dielectric constant in radial direction, 
e± = d/(d — 4:Tt(3±). Here, (3ii and /3j_ are the polarizabil- 
ities per unit area of the BN-sheet which are extracted 
from the bulk calculation 1 ^. The polarizability 7 per unit 
length of the cylinder due to an external homogeneous 
electric field perpendicular to the tube axis is2I 
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(7) 

with 9 = (R - d/2)/(R + d/2) and v = ^f^jTZ- In 
the limit R/ d — > 00, the polarizability in Eq. 10 displays 
a linear dependence on the radius: j(R) — > 7o(-R) — S, 
where 70 = +P±)R- This corresponds to the po- 
larizability without depolarization effects and coincides 
with the undamped model curve for j±_ (dotted line in 
Fig. El . 

The depolarization effects are introduced into our 
model by multiplying the undamped model curve for the 
perpendicular polarizability with the "damping" factor 
T(i?) = j(R)/'y (R). This factor depends on the cylin- 
der thickness d. The value d = 3A, which corresponds 
approximately to the full width of the charge-density of a 
BN-sheet, leads to an almost perfect agreement between 
model and ab-initio calculations. 22 

To compute Raman intensities we make the further 
assumption: 
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where da^jdu^ is constructed according to Eq. EJ We 
assume here that to first order the atomic displacement 
does not change the depolarization. For i = j = 3, i.e., 
for incoming and scattered light polarized along the tube 
axis, Tij — 1, otherwise Ty = T(R). 

In Fig. |3| we present the ab-initio and model Raman 
spectra for the (9,0), (13,0), and (16,0) zigzag BN nan- 
otubes and a (10,10) armchair tube. The latter two have 
diameters (12.8 Aand 13.8 A) in the range of experimen- 
tally produced BN tubes&S. The spectra are averaged 
over the polarization of the incoming light and scattered 
light. We first discuss the spectra of the zigzag tubes. 
The peaks below 700 cm -1 are due to low frequency 
phonon modes which are derived from the acoustic modes 
of the sheet and whose frequencies scale inversely propor- 
tional to the tube diameter (except for the E2(R) mode 
which scales with the inverse square of the diameter 
The E2(i?) mode gets quite intense with increasing tube 
diameter, but its frequency is so low that it will be hard 
to distinguish it from the strong Raleigh-scattering peak 
in experiments. The Ei(L) peak has almost vanishing 
intensity in the ab-initio spectrum and is overestimated 
in the model. The radial breathing mode (RBM) yields 
a clear peak which should be easily detectable in Raman 
measurements of BNNTs just as in the case of CNTs. 
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FIG. 3: Raman spectrum for different BN tubes: Compar- 
ison of ab-initio calculations (positive axis) with the bond- 
polarization model (inverted axis) . Symmetry assignment fol- 
lows Ref. 121 The letters R,T,L denote the character of the 
corresponding phonon oscillation: radial, transverse, or lon- 
gitudinal (see Ref. llOt) . 



Both ab-initio and model calculations yield a similar in- 
tensity for this peak. The high-frequency modes above 
700 cm -1 arc derived from the optical modes of the sheet 
and change weakly with diameter. The Ai(R) mode at 
810 cm -1 gives a small contribution which might be de- 
tectable. The intensity decreases, however, with increas- 
ing diameter. The model only yields a vanishingly small 
intensity for this peak. At 1370 cm -1 a clear signal is 
given by the Ai(T) mode which has very similar inten- 
sity both in model and ab-initio calculations. The small 
side peak at slightly lower frequency is due to the E 2 (L) 
mode. The Ei(T) peak at 1480 cm -1 is gaining intensity 
with increasing tube radius. The overall Raman spec- 
trum for a (10,10) armchair tube exhibit similar trends. 

In Fig.0]we show for the (16,0) tube the dependence of 
the intensities on the light polarizations. If both et and 
ef point along the tube axis (Fig.|3|a), only the Ai modes 
are visible and described well by the model (except the 
810 cm -1 mode). This coincides with the finding that 
for the polarizability along the tube axis, depolarization 
does not play a role24. The E modes are only visible if 
at least one of and e f has a component perpendicular 
to the tube axis (Fig. 0] b and c) . The bond polariz- 
ability model reproduces these peaks, but tends to over- 
estimate the E modes. The inclusion of depolarization 
effects is absolutely mandatory. Without depolarization, 
the model overestimates the Raman intensities for per- 
pendicular polarization by about a factor of 15. Other 
discrepancies are due to the assumption in Eq. J8J. 
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FIG. 4: Raman spectrum of a BN(16,0) tube for different 
light polarizations — » ef. Tube oriented along (OOf). 

In conclusion, we implemented the bond polarizabil- 
ity model for BN nanotubes with parameters taken 
from ab-initio calculations and under inclusion of de- 
polarization effects. Going beyond previous models for 



graphitic systems, our calculations yield different pa- 
rameters for the in-plane and out-of-plane perpendic- 
ular polarizabilities. Good agreement between model 
and ab-initio calculations of the non-resonant Raman 
spectra of BN nanotubes is obtained for light polar- 
ization along the tube axis. For perpendicular polar- 
ization, the inclusion of depolarization effects leads to 
a reasonable agreement between model and ab-initio 
spectra. The model is implemented for single-wall BN 
tubes but can be extended to multi-wall tubes if the 
strength of the depolarization effects is modeled accord- 
ingly. A similar bond-polarizability model can also be 
developed for the non-resonant Raman spectra of semi- 
conducting carbon NTs. However, due to the metallic 
behavior, a bond-polarizability model is not applicable 
to the graphene sheet. Consequently, the modeling of 
the polarizability of semiconducting tubes is very sensi- 
tive to the band-structure 25 , in particular to the band- 
gap which depends on the radius and chirality of the 
tubes. This work was supported by EU Network of Ex- 
cellence NANOQUANTA (NMP4-CT-2004-500198) and 
Spanish-MCyT. Calculations were performed at IDRIS 
and CEPBA supercomputer centers. 
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